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1. Introduction 

Tactoids are metastable nematic liquid crystalline droplets suspended in the co-existing 
isotropic fluid parent phase [TJ El lU IS1 El ■ They can adopt a variety of shapes and 
director-field configurations, depending on their size, on their elastic properties and on 
the degree of anisotropy of the surface tension [HI El HD1 [TT1 [T21 UHl E31 UHl HEl HZl [TBI - 
Theoretically, small nematic droplets tend to be quite elongated with a director field that 
is uniform, whereas very large droplets should be more or less spheroidal with a bipolar 
director field that supports two surface defects called boojums jSJ. The transformation of 
the shape and the director-field configuration of the droplets with increasing size appears 
to be continuous [§]. It should be noted that tactoids are quite different from the usual 
polymer dispersed liquid crystals or PDLCs, in which nematic droplets are embedded 
in a (solid) polymeric matrix. Indeed, in PDLCs the liquid-crystalline droplets have 
little opportunity to adjust their shape once the polymer matrix has set. Their surface 
anchoring is usually also fixed by a strongly anisotropic surface tension that is not nearly 
as small as seems to be the case in tactoids [HI HH] . 

In 1985, Williams showed theoretically that spherical bipolar droplets can also 
exhibit a parity-broken, twisted director field without the nematogens needing to be 
chiral JUj ■ Twist transitions have in fact also been predicted to occur in other geometries 
including nematic fluids confined between two concentric spheres |2D] and in so-called 
hybrid nematic layers [2TJ 122 122] • In a twisted nematic droplet the director field circles 
around the main droplet axis as shown in Figure 1. Twisted bipolar droplets have indeed 
been observed in thermotropic liquid crystals j3J|211l2H]- To our knowledge, they have 
not been observed in lyotropic liquid crystals — as we shall see, this may be connected 
with their typical elongated, spindle-like shape. 

In this paper, we extend the work of Williams by considering tactoids of finite size, 
that is, we allow for nematic droplet shapes that are not necessarily spherical, and for 
director fields that are not truly bipolar. Our aim is to estimate the minimum size for 
the twisted director field to be stable in relation to the elastic properties of the nematic 
and the interfacial tension between the co-existing nematic and isotropic phases. To 
this end, we first, in section 2, write down a phenomenological free energy consisting of 
the sum of the usual Frank energy of the elastic deformation of the director field, and 
a surface free energy of the Rapini-Papoular type [2H], and discuss the earlier work of 
Williams. 

Next, in section 3, we present our parametrization of the problem, which is based 
on one that we introduced in previous work, and that allows for (but does not impose) 
a continuous transformation from a bipolar to a uniform director field, and from a 
spherical to a highly elongated spindle-like droplet. Here, we add an extra degree of 
freedom to describe a potentially twisted director field, being the maximum surface 
twist angle, as parameter. In section 4, we first discuss truly bipolar droplets, which 
can largely be dealt with analytically, and generalize Williams' result to bipolar droplets 
of arbitrary aspect ratio. As we shall see, the larger their aspect ratio, the less likely a 
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twisted director-field configuration. We numerically investigate the more general case of 
droplets that can optimize their structure within a prescribed class of droplet shape and 
director-field geometry, in section 5. We find that the transition to a twisted structure, 
if it occurs, is sharp at a critical dimensionless volume Vq and that the droplet is twisted 
for all volumes v > v o greater than that. The critical volume vq is a complex function 
of the various elastic constants and of the anisotropic surface tension. We find that for 
lyotropic liquid crystals the twist transition must occur at droplet volumes three to six 
orders of magnitude larger than that for thermotropic liquid crystals. This presumably 
explains why twisted structures have not yet been observed in the former, whereas they 
have in the latter. Finally, in the last section, section 6, we summarize our conclusions. 

2. Free Energy 

Presuming the bulk free energy describing the thermodynamic stability of a nematic 
droplet to be independent of its macroscopic structure, the excess free energy F relevant 
to the problem in hand may be written as jHj 

F = F E + F S , (1) 

F E = \ [ d 3 r [K n (V ■ n) 2 + K 22 (n • V x n) 2 + K 33 (n x (V x n)) 2 ] , (2) 

F s = r [ rf 2 r(l + o;(q-n) 2 ), (3) 

where the volume integral F E accounts for a possible elastic deformation of the director 
field, and the surface integral F$ for the interfacial free energy. Here, n(r) denotes the 
director field at position r, q = q(r) is the normal to the droplet surface, K u , K 22 
and -KT33 are the familiar Frank elastic constants for splay, twist and bend deformations, 
r is the isotropic interfacial tension and uo > a dimensionless anchoring strength, 
penalizing a nonplanar alignment of the director field to the interface. 

In equation (J2J), we ignore the contributions of the surface elastic constants K13 
and K 24: , associated with splay-bend and saddle-splay deformations. The splay-bend 
contribution we drop for simplicity, because keeping this term would, in principle, require 
the incorporation of higher order gradients in the director field to keep the problem well 
posed |27j . We note in passing that Yokoyama recently argued on quite general grounds 
that K13 should in fact vanish identically |2Hj- As for the saddle-splay contribution, it 
merely renormalizes the splay elastic constant Ku [9J, at least for twist-free, bispherical 
director-field configurations, in which case Kn — > Kn — K 2 ^ In a twisted bipolar 
configuration this is no longer strictly accurate but we do not expect that accounting for 
the splay-bend elasticity would qualitatively affect the stability of the twisted director 
field for reasons that we discuss in more detail in section 4. It should be stressed that 
in certain geometries, such as in hybrid nematic cells with a homeotropic alignment on 
one boundary surface and planar alignment on the other, or in nematic layers with a 
mutually perpendicular planar alignment on both surfaces, the saddle splay contribution 
can induce a twist transition (2111221123], so the issue does indeed deserve some attention. 
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The optimal droplet shape and director-field configuration minimizes equation (JTJ 
at constant volume V of the drop and at constant unit length of the director |n(r)| 2 = 1. 
This highly non-trivial free boundary problem can be simplified considerably in the limit 
of V/tKhui — > oo, because then the droplet is spherical and the director field is aligned 
tangentially to the interface [H]. Tacitly presuming this limit to hold, Williams showed 
that for the twisted bipolar configuration to be favored, the elastic constants of the 
nematic would have to obey the inequality [TU] 

4tt 2 - 16, , 
733 < 2Q _ ~ 2 (1 - 722) « 2.32(1 - 722 ), (4) 

where 722 = K^/Kw and 733 = K 33 /Ku. This inequality expresses the circumstance 
that by adopting a twisted configuration, a droplet can lower its splay energy albeit 
at the cost of raising its bend and twist energy. Hence, if the bend and twist 
elastic constants are small enough compared to the splay elastic constant, a twisted 
configuration can become energetically favorable. Volovik and Lavrentovich were the 
first to observe a twisted bipolar droplet in a thermotropic nematic liquid crystal [24J. 
Available experimental data seem to indicate that equation (j3J) indeed holds |6J . 

Whilst twisted bipolar drops have been observed in thermotropic nematics, such 
is as far as we are aware not the case for lyotropic ones. This may, of course, be 
linked to them not obeying equation (J3J. Indeed, theoretical calculations show that for 
hard, rod-like particles, 722 ~ 1/3 and 733 <J 5 at least if they are sufficiently slender 
[29| l3(H 135] . These predictions agree with experimental data on aqueous dispersions 
of tobacco mosaic virus particles that to a first approximation behave like hard rods, 
with 722 ~ 0.2 — 0.4 and 733 ~ 9 — 17 [3TJ E2j- So, it seems that a nematic of rod- 
like particles interacting via purely repulsive interactions do not obey the Williams 
inequality, equation (jlj), and hence should not display a twist transition. On the other 
hand, one would expect nematic droplets in dispersions of semi-flexible particles to be 
able to support a twisted director field, because calculations show that 722 ~ 1/3 and 
733 ~ 1 |2H1 EU ES] • This is in qualitative agreement with experimental data albeit that 
the ratio of the twist and splay constant is in practice somewhat lower than predicted, 

722 » 0.04 - 0.15 mm eu mm eh eh wl- 

The question thus arises why twisted textures have not been observed in lyotropic 
nematics. To answer this question, we generalize equation (j3J) to nonspherical, spindle- 
like droplets, and allow for director-field configurations that are not necessarily truly 
bipolar. As shall become clear, both finite-size effects and non-sphericity (typical of 
lyotropic nematic tactoids) conspire against the twisted state. 

3. Parametrization 

Following the prescription of references (HI E3 El], we impose the geometry of the 
director field and that of the droplet shape. See figure 2. They are a function of 
three free parameters, being the aspect ratio of the droplet, the position of the foci of 
the bispherical director field |H1 El EDJ and the surface twist angle [H] . The director 
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field is chosen in such a way that it can transform continuously from a homogeneous 
to a bipolar field and from an untwisted to a twisted one. The droplet shape varies 
continuously from spherical to spindle like, i.e., elongated with sharp ends. For the 
given director fields and droplet shapes, discussed in more detail below, we calculate 
the total free energy, equation ( | T |) . and minimize it with respect to the free parameters. 

Inspired by the works of Williams ^Uj and Kaznacheev and co-workers ^21 E3] 
we assume the director field n to be twisted bispherical, that is, the director field lines 
connect virtual point defects [H], situated on the main-body axis of the drop and 
a distance 2R apart. (See also our previous work Each director field line lies 

on a surface of revolution of a circle section about the main axis as shown in Figure 
1. Without loss of generality, we assume the latter to coincide with the z-axis, and 
the origin to be separated a distance R from both virtual defects. As for the shape 
of the droplet, contrary to Williams' assumption, we assume it not to be necessarily 
spherical but potentially elongated. Its surface is the surface of revolution of a circle 
section around the z-axis. This shape very accurately describes the shape of actual 
tactoids [3 HH UH]- The distance between the two poles of the droplet is set equal to 
2R, where R < R, and the droplet is positioned symmetrically in the director field, 
that is, both poles are a distance R away from the origin (see Figure 2). We perform 
all the calculations in bispherical coordinates (£,7/,0), because of the symmetry of the 
problem in hand and because of the fact that a surface of constant rj is the surface of 
revolution of a circle section, see Figure 3. For further details of our parametrization 
see the Appendix A. 

We expect the twist deformation, if there is one, to be small near the line connecting 
the two poles of the droplet and to grow larger towards the surface of the droplet. The 
reason for this is that introducing twist near the center of the droplet will greatly 
increase the amount of bend deformation, whereas introducing it near the surface has 
a much smaller effect. Hence we expect the angle a between a director field line and a 
meridian on the surface of revolution on which the line resides to depend strongly on 
the ^-coordinate but much more weakly on the ^-coordinate. As a first approximation, 
we neglect the dependence on £ and assume that a = a.{rj). To keep things simple, we 
furthermore presume a to only depend on a single parameter, i.e., a (77) = ao f(v)i where 
«o sets the maximum twist angle at the surface of a drop and f(rj) is some function of 
rj. To fix the latter, we demand that ot{rj) oc 77^ if 77 J, 0, with 8 > 0. This requirement is 
needed to keep the contributions from the twist and the bend to the elastic free energy 
deformation finite, and confirms the earlier claim that near the center of the droplet the 
twist should be small. Following Williams fTU] . we use the Ansatz a{rj) = aosmr] that 
fulfills all requirements. 

Now that we have established our parametrization of the problem,we can find the 
optimal droplet shape and director-field configuration by minimizing the free energy at 
a constant volume with respect to the free parameters e, R and a . Before doing this, 
we first make the free energy dimensionless by dividing it by rV 2 ^ 3 . The dimensionless 
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free energy F then reads 

F = ^ N = F s + u;v~ 1/3 (F 11 + F 22 + F 33 ), (5) 

where v = V{tuj / Ku) 3 is a dimensionless volume, and F$ = Fs/rV 2 ^ 3 and Fa = 
Faj KuV 1 ! 3 = ^ulii/V 1 / 3 are dimensionless free energies with la = Fu/Ku] Fa is the 
contribution to Fe associated with the constant Ku. See the Appendix A. The ratios 
7n = K\\j K\\ = 1, 722 = K22/K11 and 733 = K 33 /Ku measure the deviation from the 
equal constant approximation. Because of the scaling Fa oc V 1 ^ 3 and F$ oc V 2 ^ 3 , Fa 
and Fs do not depend on the volume of the droplet, only on its shape and director field. 
The free energy equation (0) is a function of the free parameters u>, v, 722 and 733. 



4. Bipolar drops: analytical results 

As already mentioned, Williams' criterion ^Hj for the onset of the twist transition 
equation (j3J) is strictly valid for spherical bipolar droplets only, i.e., in the limit of 
infinitely large droplets jHHOj- We extend this result to droplets of finite size. Before 
turning to the more general case of quasi-bipolar director fields we first allow for non- 
spherical droplet shapes given a truly bipolar director field, which corresponds to the 
situation where uj —>■ 00. For that case we can analytically calculate a generalization of 
inequality (j3J). Note that according to theoretical estimates u ~ 1 for lyotropic systems, 
although experiments point at somewhat larger values closer to ten |%| IH| IT3j. 

Let us assume that the transition from an untwisted to a twisted director field is 
continuous, that is, the twist angle ao changes continuously from zero to a non-zero 
value, an assumption that will be supported by numerical results discussed below. Near 
the transition, we can then make a Taylor expansion of the free energy for small a>o, 



dF 

a =o da 



a + 



1 d 2 F 



2 da 2 

a=0 



at (6) 

a=0 



For symmetry reasons, the first non-constant term in cto should be of second order in cto 
so dF /da = for a = 0. At the point where the twist transition occurs, the coefficient 
of the 0(«Q)-term vanishes. This coefficient can be calculated by putting R = R (the 
director field being bipolar) and expanding (V • n) 2 , (n • V x n) 2 and (n x (V x n)) 2 
in powers of o:q, and then performing the integrations of equation (jA.ll|) . Because the 
droplet is bipolar, we have (n • q) = so the surface term Fs does not depend on a . 
We thus find for the coefficient of the 0(ag)-term 



1<9 2 F 



Ida 2 



(-4/i(6) + 4/ 2 (e) 7 22 + / 3 (e) 7 33) , (7) 



Ul/3y(e)l/3(l + e 2)2 



a=0 

where the volume V(e) is given by equation (jA.14|) . and 

f^e) = (l + e 2 ) 2 arctan 2 e-e 2 , (8) 

/2(e) = A(e) + 2e(l - e 2 ) arctane, (9) 
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/ 3 (e) = A(e)-2/ 2 ( e )+ 46 2 . (10) 

Bipolar droplets attain a twisted director-field structure, if equation (jZJ) becomes 
negative, i.e., if 

/i(e) > 722 / 2 (e) + ^733 /s(e), (11) 

because then the free energy is unstable to a perturbation away from zero twist angle. 
Equation (fTTj) is a generalization of the Williams inequality equation (@J for elongated 
nematic droplets of aspect ratio e. Setting e = 1 in equation ([TT]) reproduces equation 
(jlj), as it should. Figure 4 shows a contour plot of the reciprocal aspect ratio e as a 
function of 722 and 733 when the equality in equation (fTTj) holds. For a given value of 
e, droplets of nematic liquid crystal with a combination of elastic constants lying below 
the contour line corresponding to that value of e have a twisted structure, droplets lying 
above are untwisted. From figure 4, we conclude that when the aspect ratio increases 
(e decreases), there are fewer possible combinations of elastic constants that give a 
twisted droplet structure, so the larger the aspect ratio, the harder it is to get a twisted 
structure. Typical values of e obtainable from experiments on colloidal dispersion are 

0.1 - 0.5 0EHI31. 

As it happens, the aspect ratio of a nematic drop is not a free parameter — it sets 
itself so as to minimize the free energy. For a bipolar droplet with R = R and a(rj) = 
we find that the free energy has the form (see also [Sj 

F = 8tt(1 + e 2 )c 2 (e)c 3 (e) 2 + Siruv- 1 / 3 ^) (Jl + ^733^) c 2 (e) - 7 3 3 e arctan 2 ej , (12) 
where 



02(e) = e — (1 — e ) arctan e (13) 



C3( e )=(fy /3 (3(l + e 2 ) 2 C2 (e)-4 e 3 )- 1/3 . (14) 



and 



v 47T / 

This free energy is valid all the way to the critical conditions set by equation (fTTjl . 
To find the optimal e for a given droplet volume v and elastic constants 722 and 733, 
we simply set dF/de = 0. This produces an implicit expression for e that we do not 
reproduce here. We have only been able to find an explicit solution if -y 1 / 5 ^ -3 / 5 < 1, 
in which case e ~ (1/4)(15/87t) 1 /\; 1 / 5 cj- 3 / 5 < 1 and the drops are very elongated. For 
small e, /i(e) ~ 4e 4 /3, /2(e) ~ 2e 2 and /3(e) ~ 4e 4 , so the parity broken structure only 
occurs if 722 is very small, i.e., at most e 2 or smaller than that, unless 733 > 4/3 in 
which case the twisted state is absolutely unstable for all 722- This confirms again that 
twisted director fields are very unlikely in elongated tactoids. 

We note that since /2(e) //i(e) and /3(e) //i(e) are both monotonically decreasing 
functions of e, and since e is a monotonically increasing function of the droplet volume 
[S], that if a twist transition occurs in a bipolar droplet at v = v the droplet must be 
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twisted for all v > v but untwisted for v < v . Our numerical studies, presented in the 
next section confirm this. 

Returning to the issue whether or not the saddle-splay elastic deformation has any 
appreciable impact on the twist transition, we find that if we add its contribution 



to the total free energy, the generalized Williams inequality, equation (fTTj). becomes 



Here 724 = K-^l Kn and /4(e) = e(l — e 2 ) arctane. We immediately see that for nearly 
spherical drops we retrieve the original Williams inequality, (HJ), because /4(e) — > as 
e I 1. This implies that the splay-bend contribution is small compared to the other 
elastic contributions for nearly spherical droplets and that we can indeed neglect it. 
However, in the opposite limit where the droplet is strongly elongated this in no longer 
the case. Inserting the limit e j in equation (|16J). our generalized Williams inequality 
becomes 724 > 2722, showing that the splay- twist elastic term may indeed stabilize the 
twisted state if the splay-bend elastic constant is sufficiently large. Unfortunately, as 
far as we are aware, not even an order-of-magnitude estimate is known of 724, at least 
not for lyotropic nematics. We can, in spite of this, nonetheless make headway by 
presuming that the Nehring-Saupe equality [12] 2724 = 1 — 722 holds. (See, however, 
[2%].) Presuming it holds, our generalized Williams inequality simplifies to 722 < 1/5, 
indicating that very slender tactoids formed in dispersions of semi-flexible polymers 
could potentially undergo a parity-breaking transition induced by the splay-bend surface 
elasticity. We think that this phenomenon is unlikely to be observed in practice, because 
typically e > 0.1 00112], i-e., tactoids are in practice just not sufficiently slender. In 
addition, to become sufficiently slender they would have to be very small, and small 
tactoids tend to revert to a uniform director field as the calculations presented in the 
next section show. For these reasons, we think it justified to focus on the bulk elasticity 
of the droplets. 

5. General case: numerical results 

Now we consider the more general case, where the director field in the droplet is 
intermediate between uniform and truly bipolar. As we and others J3] recently 
argued, nematic tactoids are usually only approximately bipolar, with a director field 
that becomes increasingly uniform the smaller their size. For this more general case it is 
not possible to analytically derive a further generalization of equation (J3J). We therefore 
proceed to investigate the droplet shape and structure using numerical methods. The 
optimal values for the aspect ratio e, degree of curvature of the director field R and 
maximum twist angle «o axe determined by numerically minimizing equation (JHJ) with 
respect to these quantities at fixed values of u, v, 722 and 733, using the quasi-Newton 
algorithm E04JYF from the NAG® Fortran Library, release Mark 18. 




(15) 




(16) 
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As already mentioned, the free energy cost of introducing a twist in the director 
field near the surface of the droplet is less than that in the center. This implies that 
introducing twist near the surface of a spherical droplet is less costly in terms of free 
energy than near the surface of an elongated droplet, assuming R/R is the same for 
both droplets. We also expect that introducing twist in a bipolar droplet, where there 
is a lot of splay deformation near the poles, has a larger effect than introducing twist in 
a droplet of the same shape, but with a more homogeneous director field. 

This then leads us to the conjecture that if a transition occurs for a certain value 
v = Vq for a given set of fixed values of u, 722 and 733, the droplet must be twisted 
for all v > Vq but untwisted for v < vq. This is true for bipolar droplets as we showed 
above, but is plausible for the more general case too on account of how the aspect ratio 
e _1 depends on v. (See, again, reference One of the consequences of this conjecture 
is that if inequality (p£J) predicts a droplet to be untwisted, the director field inside the 
droplet will be untwisted at every volume. The reason is that equation (0} is formally 
valid only for v — *■ 00. 

The first case we investigate is the so called equal-constant approximation 722 = 
733 = 1, for which K\\ = K22 = ^33- This combination of elastic constants does not 
satisfy inequality so, on account of the conjecture we put forward in the previous 
paragraph, we do not expect a twisted director field in the droplet at any volume. This 
is indeed what we find from the numerical minimization: for uj = 0.01, 1 and 100, cto = 
for all v tested in the range from 10 -18 to 10 18 . Next, we study a combination of elastic 
constants valid for rod-like colloids. Putting 722 = 1/3 in equation (jlj), we see that for 
a spherical bipolar droplet to be twisted, we must have 733 < 1.55, which is much lower 
than the typical value for rod-like colloids of 733 ^ 5. (See section 2.) So, again from 
the conjecture from the previous paragraph, we do not expect a twisted director field 
at any droplet volume. This is confirmed for 733 = 5 and to = 0.01, 1 and 100 by a 
numerical minimization for dimensionless volumes from 10~ 18 to 10 18 . 

Let us now keep the value of 722 fixed at 1/3, but take for 733 a value of 1. Inequality 
now predicts a twisted director field to be stable in the limit of an infinite volume 
implying that it might also be stable for some finite volume. That this is indeed the 
case, is shown in Figure 5, where we show the twist angle «o as a function of the 
dimensionless volume v for uj = 0.01, 1 and 100. We see that the transition is sharp and 
continuous, which was our assumption in deriving the results of section 3. This means 
that the droplet is absolutely untwisted until the volume reaches the critical size v = vq, 
from whereon the twist gradually increases with increasing v. For u = 100, the twist 
transition is not shown, but according to our calculations it occurs at 10 logt>o — 10.60. 
This can in fact be checked numerically via a different route. The aspect ratio at the 
transition can be calculated using the generalized Williams condition equation 
Because the director field at the transition is nearly bipolar (see figure 6 and the next 
paragraph), the aspect ratio should also minimize equation (|12p. which then determines 
vq. We find that 10 logfo = 10.61 in good agreement with our previous estimate. 

In figure 6 we have plotted a diagram of states for 722 = 1/3 and 733 = 1. The 
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twist transition is the only transition in the diagram that is sharp, i.e., on one side of 
the line demarcating the transition, the droplets are untwisted whereas they are twisted 
on the other. The line demarcating the transition from a homogeneous to a bispherical 
director field is drawn at the point where R = 2R, and one showing the transition in the 
aspect ratio from spheroidal to elongated is set at e = 1/2. Both are smooth crossovers, 
discussed in more detail in reference [UJ. This means that a droplet, e.g., in that part 
of the phase diagram indicated by "spheroidal twisted bipolar" is not exactly spherical 
nor exactly bipolar, but nearly so albeit that the director field is indeed twisted. 

Figure 6 shows that, at the twist transition, v is independent of uj when uj is small, 
whereas for large values of u this is not so with v ~ uj 3 . The first result can easily be 
understood by noting that ui describes the degree of anisotropy of the anisotropic surface 
tension. If uj is small, this anisotropy is small and it will not influence the droplet shape 
nor the director-field configuration. The second observation follows from the fact that 
when the droplet is bipolar, the aspect ratio must be constant at the twist transition. 
This can be seen from equation pijl . which is independent of uj, implicating that the 
aspect ratio is constant as a function of uj if 722 and 733 are constant. We have shown 
in previous work [Oj that if uj 3> 1, the aspect ratio e of a quasi-bipolar droplet scales 
as v -1 / 5 c<j 3 / 5 , from which one can conclude that v ~ uj 3 if e is constant and uj 3> 1. (See 
also equation (JUJ).) 

Next, we take 722 = 0.05 and 733 = 1, which are representative values of the elastic 
constants for semi-flexible particles (see section 2). In figures 7 and 8, we show the twist 
angle «o as a function of the dimensionless volume v for anchoring strengths uj = 0.01, 
1 and 100, and the corresponding diagram of states. The behavior of the droplets is 
similar as for the previous set of parameter values, except that the twist transition now 
occurs at a much lower value of v and that the diagram of states is somewhat more 
complex. This is to be expected, considering that the twist elastic constant is much 
smaller. The twist transition for a strong anchoring with uj = 100 is not shown, but it 
occurs at 10 logt'o = 7.71. Using the same method as described earlier in this section 
we can check this result with the theory of section 4. We find 10 logt>o = 7.79. The 
discrepancy between the two results is probably caused by the fact that the droplet is 
not quite bipolar yet (see figure 8). 

As already mentioned, we expect uj ~ 1 to be a reasonable estimate for the 
dimensionless anchoring strength. (See again, however, [SJ|UJ[T3].) Because the elastic 
constants and their ratios differ from system to system, we show, in Figure 9, the 
dimensionless volume v at which the twist transition occurs for various values of 722 
and 733 at a fixed anchoring strength uj = 1. As shown in the figure, the larger 733, the 
larger the volume of a tactoid needs to be to be able to support a twisted director field. 
The reason is that a twist deformation goes hand in hand with a bend deformation that 
is penalized more the larger the bend elastic constant K33. Also, the larger 722, the 
larger the elastic free energy penalty associated with a twist deformation, the smaller 
733 must be for the twisted director field to be favorable. The shift to higher droplet 
volumes is, as discussed, linked with the greater ease of a twist deformation the larger 
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the drop is because large drops are more spherical than small ones. It is useful to point 
out that provided the values of the elastic constants and the surface tension are known, 
figure 9 can be used to estimate the minimum droplet volume at which a transition to 
a twisted structure should be expected to occur. 

We again turn to the question why twisted nematic droplets have so far evaded 
observation in lyotropic liquid crystals whereas they have been observed in thermotropic 
ones [3J |21]. We have seen that we should not expect twisted nematic droplets in 
nematics of hard rods because 722 and 733 are too large. For worm-like, semi-flexible 
particles such as fd virus though, 722 and 733 possibly have the right values but the 
volume of the droplets has to be sufficiently large for parity breaking to occur. 

In fact, we can use figure 9 to estimate how large. First, we know that v oc Vr 3 . 
Since the elastic constants of a thermotropic and a lyotropic crystal are of the same 
order of magnitude, we expect the twist transition to occur at about the same value of 
v. However, r is typically one or two orders of magnitude larger in a thermotropic liquid 
crystal than in a lyotropic one [TU ESI EH ESI ESI EZl EH! EOl EOJ - This, then, means that 
the volume of a droplet of lyotropic liquid crystal at the twist transition must be three 
to six orders of magnitude larger than that of a thermotropic liquid crystal. 

To make our estimate more concrete, we set K 33 /K n = 1, K 2 2/K n = 0.05 and 
u = 1, and read off from figure 9 that the twist transition occurs if Vq = V(tuj/Kh) 3 = 
10 1 ' 35 . Inserting the typical values Ku = 10~ n N and r = 10~ 6 N/m, this corresponds 
to a droplet volume of V ~ 10 4 /im 3 , with a linear size of about 30 /im. However, if uj 
is closer to 10, as might well be the case 01131111!, this value increases to > 300 fim. A 
similar enhancement occurs if r is larger by a factor of ten. 

6. Conclusion 

By adopting a twisted director field structure, a nematic tactoid reduces its splay energy 
albeit at the expense of increasing its bend and twist energy. Hence, if the bend and 
twist elastic constants are small enough in comparison to the splay elastic constant, 
this twisted (parity-broken) structure is energetically more favorable. For a spherical 
bipolar droplet, Williams [TU] derived a criterion, equation (jlj), the elastic constants 
have to obey in order to be able to support a twisted director field. In this work we 
generalized this criterion to bipolar droplets of arbitrary aspect ratio, equation (JTTJ), 
showing that the more elongated the drop the less likely they are to have a twisted 
director field. 

We find the transition to a twisted structure to be sharp, in contrast to the 
transformation from a uniform director field to a bipolar one. When the transition 
occurs at a certain dimensionless droplet volume v = vo, the droplet is twisted for all 
v > v but untwisted for v < v . One of the consequences is that droplets that do not 
have a twisted structure at infinite volume (the limit in which Williams' result is valid), 
cannot have a twisted structure for any finite volume either. Not entirely surprisingly, v 
increases with increasing 733 = K 33 /Kn and decreases with increasing 722 = ^22/^11 ■ 
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Nematics of rod-like colloidal particles appear to have elastic constants that do 
not obey the Williams inequality, so we do not expect to find droplets with a twisted 
structure in this sort of system. This is not so for dispersions of semi-flexible colloidal 
particles (or polymers), however, but the tactoids have to be large enough for the twisted 
director-field configurations to become stable. This is in part due to them becoming 
less elongated with increasing size. 
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Appendix A 



The bispherical coordinates are related to the Cartesian coordinates x = (x,y,z), 
through the relations 

x = RZ~ X sin?7sin£cos</>, 

y = — RZ^ 1 sin 7/ sin £ sin 0, (A.l) 
z = i?Z _1 cos£, 

where Z — 1 + sin £ cos?] and < £ < 7r, < r] < ir and < <fi < 2tt. The Cartesian 
unit vectors e x , e y and e z are related to the bispherical unit vectors eg, e v and e^ by 

e x = Z~ x sin 77 cos £ cos eg + Z -1 (sin£ + cos 77) cos</> e^ — sin0 e^, 
e y = — Z^ 1 sin 77 cos£sin</> eg — Z _1 (sin^ + cos?]) sin0 e,, — cos0 e^, (A. 2) 

e z = — Z _1 (sin^ + cos?]) eg + Z^ 1 sin?7cos£ e^. 

Finally, /ig = |<9x/<9£| = RZ' 1 , h v = \dx./dri\ = RZ" 1 sin£ and = \dx/d(j)\ = 
RZ~ l sin £ sin 77 are the metric components of the bispherical coordinate system. 
The director field n is given by 



where 
and 



do 



cos a (77) eg 



d 1 = (R 2 — R 2 )Z sin r\ cos£, 



cos airf)^ + sina(7/)e (? 



(A.3) 

(A.4) 
(A.5) 



d 2 = (R 2 -R 2 )Z cost] + {R 2 + R 2 )Zsin£. 

This choice of director field leads to the following expression for the various components 
of the elastic free energy, equation Q, 

Vn = — , = = (4 cos£ cos a + dia v sin a), (A. 6) 



^/dj+d] 
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Z - 1 + d 2 a, , (A.7) 



x (V x n)) 5 = cot £ sin 2 a + — 



Zdi cos a 
sin£(<i 2 + dl) 



(2 sin 2 £ sin r\ cos a — c^a:^ sin a) 



(n x (V x n))„ 
and 

(nx(Vxn)) (f = 



sin 2 a (sin £ + cos 77) d\Za ri sin2a 2d 2 Z sin £ sin 77 cos 2 a 



+ 



+ 



sin £ sin 77 ' 2 sin£(<i 2 + d|) ' d\ + d \ 

cos a f d\ sin a (sin £ + cos?]) diZa^cosa 



(A.8) 
(A.9) 



+ 



— d 2 cot £ sin a , (A. 10) 



sin £ sin 77 sin £ 

where = da{rj)/dr] and (n x (V x n))j denotes the component of n x (V x n) in the 
z-direction, where i G {£,77,0}. 

Within the bipolar coordinate system, the surface of a droplet with aspect ratio 
e -1 is given by 77 = r] = 2arctane, < £ < tt and < < 2tt. So, using equations 
(jA.6j) - ([A.lOp . we can calculate the splay, twist and bend contributions F u = i^n/n, 
F22 = K22I22 and F 33 = -^33/33 to the Frank elastic energy, F E = F u + F 22 + F 33 of 
equation (J2J), where 

In = \ f v dh (V • n) 2 = \ J* d V P Q ° di d^h^ (V ■ n) 2 , 

I22 = lJ v d 3 r(n-V x n) 2 = \ £ drj J™ d£ d^h^h^ (n • V x n) 2 , 

J 33 = \ f v dh (n x (V x n)) 2 = § ft d V JZ° d£ fyhthnht (n x (V x n)) 2 
The surface free energy becomes 

F s = t d 2 r(l + lu(ii ■ q) 2 ) = r d£ d<f>{\ + uj(n ■ q) 2 )/^!^, 
Js Jo Jo 

with, for our choice of director field and droplet surface, 

dj 



(All) 



(A.12) 



cos a. 



dj + dl 

For the volume V and the surface area S of the droplet we find 



(A.13) 



V(e) 
and 



2tt 



7T 



^7 / ^£ / dcph^hrjh^ = —R 3 



1 + e 



2\ 2 



1 



1-e 2 



arctan e — 4 



(A.14) 



5(e) 



2- 



d£ / d^hthflr,^ = 2nR 



2 l + e 2 



1-e 2 
1 arctan e 



.(A.15) 
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Figure 1. Sketch of a twisted director-field line in a spindle-like nematic droplet 
or factoid. The twist is largest near the surface of the drop. Close to the center line of 
the drop the twist is presumably negligible. 

Figure 2. Droplet surface (thick lines) and defining surfaces for the director field 
lines (thin lines). The system is rotationally symmetric around the line connecting the 
two (virtual) defects given by the crossing or focal points of the director field lines. R 
is half the length of the major axis of the droplet. R is half the distance between the 
virtual defects. 

Figure 3. The bispherical coordinate system. The unit vector e^ points out of 
the plane on the right hand side of the line indicated by r\ = 0. The picture should be 
rotated along the line rj = to get the full three dimensional coordinate system. See 
the main text for an explanation of the symbols. 

Figure 4. Contour plot of the inverse of the aspect ratio e of a factoid exactly 
at the twist transition as a function of the ratios 733 = K^/Ku and 722 = K22/K11 
for bipolar droplets where if 33, if 22 and ifn are the usual bend, twist and splay elastic 
constants. The numbers next to the contour lines denote the values of e. A bipolar 
droplet has a twisted structure if and only if the point denoting the combination of 
the ratios of elastic constants 722 and 733 of the nematic lies below the contour line 
corresponding to the aspect ratio of the droplet. 

Figure 5. The maximum (surface) twist angle a as a function of the dimensionless 
volume v in units of it radians for the ratios of the elastic constants 733 = K 33 /K u = 1 
and 722 = if 22 /if 11 = 1/3 at three different values of the anchoring strength to. For 
uj = 100, the twist transition occurs at 10 logt> = 10.60 (not shown). 

Figure 6. The diagram of states for 733 = K^/Ku = 1 and 722 = K22/ Ku = 1/3. 
The diagram shows the shape of the droplet as a function of the dimensionless volume v 
and the anchoring strength u. Indicated are the transition from a uniform director field 
to a bipolar one (vertical line), the transition from a spheroidal to an elongated shape 
(horizontal line) and the transition from a bipolar to a twisted director field (line in the 
lower right corner). Only the latter transition is sharp, the transition from uniform to 
bipolar is (arbitrarily) taken at the boojum position R/R = 2 and the transition from 
spheroidal to bipolar at an aspect ratio of 2. 

Figure 7. The twist angle a as a function of the dimensionless volume v in units 
of 7r radians for 733 = K^/Ku = 1 and 722 = if 22 /if 11 — 0.05 at three different values 
of the anchoring strength tv. For uj = 100, the twist transition occurs at 10 log-y = 7.71 
(not shown). 

Figure 8. As figure 6, but now with 733 = if 33 /ifn = 1 and 722 = if 22 /if 11 = 0.05. 
The various director-field patterns and droplet shapes are again indicated. Notice the 
increased number of regimes. 
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Figure 9. The value of v at which the twist transition occurs as a function of 
733 for various values of 722 = -^22/^11 and with uj — 1. For given 722, the droplets 
above the line have an untwisted structure, whereas the droplets below it have a twisted 
structure. By adopting a twisted structure, a droplet decreases its splay energy, at the 
same time increasing its twist and bend energy. For that reason, the twisted structure 
occurs only if 733 = K 33 /K n is small enough, even when 722 = 0. 
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